Numerical evidences of spin-1/2 chain approaching spin-1 chain 
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Abstract 

In this article, we study the one dimensional Heisenberg spin-1/2 alternating bond chain in which 
the nearest neighbor exchange couplings are ferromagnetic (FM) and antiferromagnetic (AF) al- 
ternatively. By using exact diagonalization and density matrix renormalization groups (DMRG) 
method, we discuss how the system approaches to the AF uniform spin-1 chain under certain con- 
dition. When the ratio of AF to FM coupling strength a {a = Jaf/Jf) is very small, the physical 
quantities of the alternating bond chain such as the spin-spin correlation, the string correlation 
function and the spin density coincide with that of the AF uniform spin-1 chain. The edge state 
problem is discussed in the present model with small a limit. In addition, the Haldane gap of the 
AF uniform spin-1 chain is 4-times of the gap of the system considered. 
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Haldane argued that the elementary excitation spectrum in the one dimensional (1-D) an- 
tiferromagnetic Heisenberg (AFH) spin system with integer spin is gapful (massive) whereas 
for the system with half-integer spin is gapless (massless) . P] Furthermore, the spin-spin cor- 
relation in the ground state for 1-D AFH spin-1 systems should have the exponential decay 
behavior due to the gapful excitation energy. By contrary, in 1-D AFH spin-1/2 systems, the 
spin-spin correlation function exhibits power law decay behavior due to the gapless excitation 
energy. There are many numerical studies to investigate this conjecture. On the 

other hand, the spin gap has been observed in the compound NENP[Ni(C2H8N2)2N02C104] 
through the susceptibility measurement by inelastic neutron scattering experiment. [21] 

Among different studies of the Haldane problem in AFH spin-1 chain, the interested one 
is given by AKLT model. ^ This model gave a simple picture for the ground state of 1-D 
AFH spin-1 system. Each S=l spin was viewed as a triplet state of two spins with S = 1/2. 
According to this picture, only two end spins of the open AKLT chain are left and form two 
"free" S* = 1/2 spins which are called the edge state and all the other = 1/2 spins form 
RVB type spin singlet state between nearest neighbor spins. In the thermodynamic limit, 



d degeneracy due to the two S =1/2 spins. This picture was 
and discussed theoretically. J]|ll| For an open spin-1 chain 



the spin-1 chain is fourfo 
observed experimentally 9] 
the Haldane gap is the difference between the lowest eigen-energy in subspace Sf^^ = (or 
^tot ~ 1) ^"^^ ^tot ~ 2- Although spin-1 chain is a disorder system, it still has a hidden 
order, the corresponding order parameter proposed by M.den Nijs and K. Rommelse [l^ 
is represented by the string correlation defined in the later discussion. Since after, the 
numerical studies have been performed for 1-D spin-1 Heisenberg model and have proven 



that the systems really posses the hidden order. 

In this article the Haldane problem is studied from another reversed approach. We start 
from a 1-D AFH spin chain with S = 1/2, the exchange couplings between nearest neighbor 
spins take Ji and J2 alternatively as shown in Fig.l. The Hamiltonian is 

H = Ji ^ Si - S + J2 ^ S j ' S j+i (1) 

<i£odd> <jrget)en> 

This model which is called alternating bond chain has been studied in detail. |13] In general, 
the excitation spectrum of this model is gapful except the case of Ji = J2, i-e. it reduced to 
the uniform spin-1/2 chain. We may ask the questions that, what is the relation between the 
Haldane gap and the gap of model (1), and how can we understand the edge state of AKLT 



model in a reasonable way. For definiteness, we are interested in the special case, where 
Ji < (ferromagnetic coupling FM) and J2 > (antiferromagnetic coupling AF). We define 
a parameter a = \J2/Ji\ and explore how the magnetic behavior of the system evolutes 
with the parameter a. We use two numerical techniques including the exact diagonalization 
method up to the system N=32 (N is the number of sites.) with periodic boundary condition 
and DMRG method^lG^ to even more large system (N=360) with open boundary condition. 

At first, in the case of large a limit, the two nearest neighbor spins connected with AF 
bond will form an AF singlet dimer and the 1-D chain reduces to N/2 S=0 dimers connected 
with FM bond. As well known that the first excited state is the state in which one of the N 
S=0 dimers is excited to become a triplet, i.e. the excitation energy is A = J2. In opposite 
case of small a limit, the two nearest neighbor spins connected with a FM bond will form 
a FM triplet dimer and the AF coupling exists between triplet dimers. Thus we can image 
that the system will behave like an AFH uniform spin-1 chain. In the following discussion 
we will argue this imagination through some physical properties calculation for Hamiltonian 
(1) in the small a cases and compare them with corresponding properties of the 1-D AFH 
spin-1 chain. 

We find that all the properties for the two systems coincide with each other in the small a 
limit. In Fig. 2, the exact diagonalization result of the model (1) for the spin-spin correlation 
on 32 sites is shown by solid circles. The plateaus in the spin-spin correlation occur at the 
pairs of site. Within each pair the spins couple each other with the FM interaction. This 
means that in small a limit each pair of spins {S = 1/2) corresponds to a FM triplet dimer 
as mentioned above. The spin-spin correlation of the AF uniform spin-1 chain on 16 sites is 
also shown in the same figure for comparison. (The data of spin-1 chain should be reduced 
by a factor 4 in comparison with the data of model (1).[5]) It can be seen that when a = 0.01 
the locus of plateaus of model (1) approaches to that of the AF uniform spin-1 chain. 

Beside the spin-spin correlation, the string correlation of the alternating bond chain also 
shows Haldane-like behavior. We define the string correlation of model (1) as follows 

^IhU) = ((.St + S^,) eM^^Y.iS^2k-i + S^,MS2,-i + 52,)|) (2) 

where (..) means taking the average over the ground state of the model considered. This 
is a simple generalization of the definition of the string correlation of integer spin S 
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chain [3] 

C/5(j) = (5iexp[z^5^5|]5;) (3) 



fc=i 

In considering the two spins with FM couphng in a S" = 1/2 alternating bond chain will form 
a triplet dimer in small a limit, which is corresponding to a spin-1 chain. The numerical 
results of the string correlation for two models are shown in Fig. 3. We can see that upon 
decreasing a, the string correlation of the alternating bond chain is closer and closer to that 
of the AF uniform spin-1 chain. 

In addition to the spin-spin and the string correlations, the average of z-component of 
spin at each site I^Sf) of the alternating bond chain will also show the "edge state" with 
average spin at the "end sites" proposed by AKLT model. In the following discussion we will 
understand what do the "edge state" and "end sites" mean in our approach: From 5 = 1/2 
model (1) to an uniform spin-1 chain. We now calculate the (5*^^) of model (1) in subspace 
Sfot = 1 by using DMRG method and consider open boundary condition for both systems. 
In Fig. 4 (a) the results are shown for the system with N=60 lattice sites and a = 0.01. It 
is clear that the two spins with the FM coupling are pairing to triplet dimers, while the 
coupling between nearest neighbor dimers is AF. The average value (S"?) for the two spins 
in each dimer are nearly the same, e.g. each spin of the edge dimer gives {Si 2) — 0.27. 
Upon entering into the interior of the chain, the (S'f ) decreases. At the two centre dimers, 
the (Sf) are zero. Increasing the size of the chain to N=120, we can see from Fig.4(b) that: 

i) the average value of (Sf ) in each dimer is nearly equal to the corresponding one in the 
case of N=60, and takes nearly the same value irrespective to a = 0.01 or a = 0.001. The 
loci of dimers are well-coincident with the (S'f)/2 of the uniform spin-1 system, e.g. for 
a = 0.001, the {Sf + Sf) ~ 0.531 for the edge dimer of model (1), it's considerably close 
to the value (0.532) of the end spin in the AFH spin-1 chain calculated by S.R.White, jiol 

ii) more dimers of spins located in the middle part of the chain with (Sf ) ~ 0. From these 
results, we conclude that the "edge state" in the AFH spin-1 system is actually the state of 
edge dimer in model (1). 

In addition, from Fig.4(c), we can see that as the N increases further, such as N=240 
and 360, the two side regions where the (Sf) 7^ change little because they come from 
the "surface" effect. So, we may image that, in the limit of 00, Nq/N 1. {Nq is 

the number of sites in the middle part.) The effect is due to the spin disturbance which 



is happened at the surface (surface mode) and the penetration depth ^ of the disturbance 
is independent of the system size, but it will depend on a. The decay behavior of (5*^) 
away from the end point of the chain in model (1) as shown in Fig. 4(c) can be fitted by the 
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exponential form \ {Sf)\ = |(5'Q)|e^«. In the small a region, decreases upon increasing a 
and the relation between ^ and a is approximately linear. Through numerical calculations, 
we could conclude that ^~12asA^— >oo and a ^ 0. The penetration depth ^ is 
approximately the double of the decay length of the AFH spin-1 chain and this is 
because two spins connected with a FM coupling of model (1) correspond to a S=l spin of 
the AFH spin-1 chain. These conclusions support AKLT model but it should be noted that: 
i) the edge state in which the average value of end spin equal to nearly S/2 is valid for both 
the AFH uniform spin-1 system and S = 1/2 system of model (1) as state above, i.e. it is 
not necessary explained by AKLT model in general, ii) The side region, originating from 
"surface" effect, is finite and independent of the system size. 

At last, we turn to study the excitation energy A1/2 of model (1) in small a limit, A1/2 is 
defined as A1/2 = Ei— Eq, where Eq and Ei are the energies of ground state and the lowest 
excite state of model (1) with periodic boundary condition respectively. For open boundary 
condition, one must define the gap A1/2 (Ai) of the alternating bond chain (the uniform 
spin-1 chain) as the difference between the lowest-lying states with St = 2 and St = (or, 
St = 2 and St = l)-0 3 ^ matter of fact, the two gaps of model (1): A1/2 = E2 — Eq 
and Ai/2 = E2 — Ei are equal. The A for both the 5 = 1/2 model (1) with Ji = 1, 
J2 = <y = 0.002 and the AFH spin-1 system with the same J2 (i.e. The Hamiltonian of the 
spin-1 chain is H = Si ■ S* j+i, here J = J2) are shown in Fig.5. In the thermodynamic 
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limit (A^ — » 00), the excitations of both systems through above definition are gapful. (For 
the spin-1 chain, Ai = AnaidaneJ, ^Haidane ~ 0.41. Pjsi) This is the so called Haldane gap 
and, however, the gap of the uniform spin-1 system is nearly 4-times of the gap of model 
(1). If we recalculate Ai for the AF uniform spin-1 system with J = J2 = a/4 and its locus 
is denoted by stars, then we find that they nearly complete coincide with the corresponding 
locus of model (1). The appearance of factor 4 is due to the following reason. The ground 
state energy of model (1) is 

= (0|i/|0) = Ji (S^■^^+l)0 + J2 J] ( 5, " 5,+i)o (4) 

<iGodd> <j&even> 

where |0) is the ground state. 



From Fig. 6(a) we can see that the mean value of local bond strength {Si ■ SiJ^i) 3=1/2 of all 
pairs of spins coupled with a FM interaction equal to 0.25 and contributes an energy 0.25Ji. 
In same, the mean value of all pairs of spins coupled with a AF interaction is equal to —0.35 
and contributes an energy of 0.35 J2- In small a limit, all the pairs with FM coupling form the 
triplet dimers, and between the dimers, there exist weak AF interactions. Now let us consider 
the lowest excite state. It is obviously that if we change the spin state of a triplet dimer as a 
whole, the energy cost is about 4 x 0.35 J2- This is much lower than the excitation energy in 
which only single one spin state is changed. This is about 2(0.25Ji + O.35J2). So the gap of 
model (1) should be approximately written as A1/2 — 4 x 0.35J2. In similar consideration, 
the excitation energy of the AF uniform spin-1 system is Ai ~ 4:J{Si ■ Si^i)s=i = 4 x I.4OJ2 
(see Fig. 6(b) Therefore we have got the relation Ai = 4Ai/2 or equivalently the J2 in 
the spin-1 system effectively equals to 4 times of J2 in model (1). In addition, like the AFH 
spin-1 systems, the excitation of model (1) between the lowest-lying St = 1 and St = 
states is gapless due to the four-fold degenerate ground state generated by two free spins of 
ends of systems. 0] jj] [3| 

In summary, we presented numerical results of model (1) and compare with the AFH 
uniform spin-1 chain. From the spin-spin correlation and the string correlation, the system 
with small a has Haldane-like behavior. The plateaus of the spin-spin correlation and the 
spin density of model (1) with small a display that spins coupled with a FM interaction 
form a FM triplet dimer. The edge state of model (1) is also observed and (5*^ + 5*2) of 
model (1) is nearly equal to (Sf) of the AFH spin-1 chain. Furthermore, the penetration 
depth ^ of model (1) is the double of the decay length of the AFH spin-1 chain and then 
the lattice constant of AFH spin-1 chain is equivalent to the double of that of model (1). 
In addition, we find that the gap of model (1) with small a is not only finite but it's also 
a quarter of the Haldane gap of the AFH uniform spin-1 chain. This is because the energy 
cost for changing the state of the triplet dimer as a whole is less than that for changing only 
one spin of the triplet dimer and the local bond strength of the AFH uniform spin-1 chain 
is 4-times of that of model (1). 
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Figure 1: The structure of an alternating bond chain. 

Figure 2: The spin-spin correlation of the ground state of the alternating bond chain 
(N=32) with periodic boundary condition. It shows many plateaus. In each plateau, the 
two spins have the same spin correlations. 

Figure 3: The string correlation of the alternating bond chain (N=32) with periodic 
boundary condition. 

Figure 4: The spin density of an alternating bond chain in subspace S^^^ = 1 in two 
cases: (a) a = 0.01, = 60, (b) a = 0.01 and a = 0.001, = 120. The data (5/)/2 
corresponds to the AFH spin-1 chain with N=60. (c) The comparison between several cases 
with different numbers of sizes, each locus corresponds to a half of different N. 

Figure 5: A1/2 is the gap of model (1) with open (or periodic) boundary condition and 
a = 0.002; Ai is the Haldane gap of the AF uniform spin-1 chain with J = J2 = 0.002 (or 
0.0005). 

Figure 6: (a) The comparison of the local bond strength L{i) of the ground state between 
model (1) of 120 sites with a = 0.001 and the AF uniform spin-1 chain of 60 sites. L{i) are 
{Si ■ Si+i)s=i/2 and {Si/2 ■ Si/2+1) s=i/^ for model (1) and the AFH spin-1 chain, respectively, 
(b) The local bond strength {Si ■ Si+i)s=i of the ground state of the AF uniform spin-1 chain 
with N=60. 
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Fig.3 H.H. Hung & CD. Gong 



0.3 - 



alternating bond chain, Open BC 



0.2 - 



a^O.01, N^60, S^^^^l 



0.1 - 



A 

N 

CO 
V 



0.0 - 



■0.1 - 



■0.2 







10 



20 



30 



40 



50 



60 



Fig.4 (a) H.H. Hung & CD. Gong 
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Fig.4 (b) H.H. Hung & CD. Gong 
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Fig.4 (c) H.H. Hung & CD. Gong 
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